Long-term Phase-coherent X-ray Timing of PSR B0540— 69 



Margaret A. Livingstone ^, Victoria M. Kaspi, Fotis P. Gavriil 

Department of Physics, Rutherford Physics Building, McGill University, 3600 University 

Street, Montreal, Quebec, H3A 2T8, Canada 



ABSTRACT 

We present a new phase- coherent timing analysis for the young, energetic 
pulsar PSR B0540— 69 using 7.6 yr of data from the Rossi X-Ray Timing Explorer. 
We measure the braking index, n = 2.140 ± 0.009, and discuss our measurement 
in the context of other discordant values reported in the literature. We present 
an improved source position from the phase- coherent timing of the pulsar, to 
our knowledge, the first of its kind from X-ray pulsar timing. In addition, we 
detect evidence for a glitch which has been previously reported but later disputed. 
The glitch occured at MJD 51335 ± 12 with Aiy/u = (1.4 ± 0.2) x 10"^ and 
Au/u = (1.33 ±0.02) X 10~^. We calculate that the glitch activity parameter for 
PSR B0540— 69 is two orders of magnitude smaller than that of the Crab pulsar 
which has otherwise very similar properties. This suggests that neutron stars of 
similar apparent ages, rotation properties and inferred dipolar B fields can have 
significantly different internal properties. 

Subject headings: pulsars: general — pulsars: individual (PSR B0540— 69) — X- 
rays: stars 



Introduction 



Young, rotation-powered pulsars spin down sufficiently rapidly to provide interesting 
tests of the simple model of pulsar spin down, given by 

z> = -Ku"", (1) 

where u is the spin frequency, z> is the frequency derivative, and K is a constant related to 
the magnetic field, moment of inertia, and angle between the spin and magnetic axes. The 
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braking index, n, is equal to 3 for magnetic dipole radiation in a vacuum. A measurement of 
the second frequency derivative, i), allows for the calculation of n. Measurements of n give 
crucial insight into the electrodynamics of pulsars. Only five such measurements have been 
made (Lyne et al. 1993, 1996; Dcctcr et al. 1999; Camilo et al. 2000; Livingstone et al. 2005). 
Interestingly, all of these values arc less than 3 and span a relatively wide range: n = 1.4±0.2 
to n = 2.91 ±0.05. The physical implications of n < 3 are still debated; possible explanations 
include the removal of angular momentum by a pulsar wind (Manchester & Peterson 1989) 
and a time-dependent magnetic field (Blandford & Romani 1988). 

Some pulsars, especially the youngest ones, exhibit significant deviations from the stan- 
dard spin-down model. Two types of deviations arc known, timing noise - a low-frequency 
random process of unknown physical origin superposed on the deterministic spin-down of 
the pulsar - and glitches, sudden increases in rotation rate likely caused by catastrophic 
unpinning of superfluid vortices (e.g. Lyne & Smith 1990). 

One of the objects for which n has been reported, ~ 1500 yr old PSR B0540— 69, has had 
several contradictory values quoted in the literature. The measured values were obtained 
from timing observations made using a variety of instruments operating from radio frequen- 
cies into the X-ray energy range and have produced values that range from n — 1.81 ± 0.07 
to 2.74 ± 0.10 (Zhang et al. 2001; Ogclman & Hasinger 1990). The differences among the 
measured values clearly signal contamination by one or more of timing noise, unidentified 
glitches, or even pulse counting errors. 

In this paper, we re-examine 4.6 yr of data from the Rossi X-ray Timing Explorer 
(RXTE) reported on by Zhang et al. (2001) and Cusumano ct al. (2003) who arrived at 
contradictory conclusions. We also extend the data set to a total of 7.6 yr. We present a 
measurement of n for PSR B0540— 69 with a partially phase-coherent analysis, and confirm 
it by performing a fully phase coherent analysis. We show that there is indeed evidence for a 
glitch near MJD 51335 ± 12 as reported by Zhang ct al. (2001), but disputed by Cusumano 
et al. (2003). However, we show that the correct value of n is 2.140 ± 0.009, within 2a of 
that reported by Cusumano et al. (2003). 

2. Observations 

Observations of PSR B0540— 69 were obtained using the Proportional Counter Array 
(PCA; Jahoda et al. 1996) on board RXTE. The PCA consists of an array of five col- 
limated xenon/methane multi-anode proportional counter units (PCUs) operating in the 
2-60 keV range, with a total effective area of approximately 6500 cm^ and a field of 
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view of ~ 1° FWHM. We used 7.6 yr of public archival RXTE observations collected in 
"GoodXenon" mode, which records the arrival time (with l-//s resolution) and energy (256 
channel resolution) of every unrejected event. We used all layers of each PCU in the 2-18 keV 
range. A similar energy range was used by both Zhang et al. (2001) and Cusumano et al. 
(2003), as it maximizes the signal-to-noise ratio for this source. Because PSR B0540— 69 was 
most often not the primary target of RXTE in these observations, integration times range 
from ~1 to >10ks, resulting in a variety of signal-to-noise ratios for individual pulse profiles. 

The observations were reduced using standard FITS tools as well as specialized soft- 
ware developed independently. Data from the different PCUs were merged and binned at 
1/1024 ms resolution. The data were then reduced to barycentric dynamical time (TDB) at 
the solar system barycenter using, at first, a position determined by Chandra observations 
(Kaaret et al. 2001) and the JPL DE200 solar system ephemeris. Once a timing solution 
was determined, we were able to fit for the position of the pulsar, as explained in §3.3. We 
then re-baryccntcrcd the data at the best-fit timing position. Each time series was folded 
with 32 phase bins using the ephemeris obtained by Deeter et al. (1999), as this ephemeris 
produced the highest quality pulse profiles. Resulting profiles were cross-correlated in the 
Fourier domain with a high signal-to-noise ratio template created by adding phase-aligned 
profiles from all observations. We implemented a Fourier domain filter by using the first 
6 harmonics in the cross-correlation, appropriate for this broad pulse profile. The cross- 
correlation produces an average time of arrival (TOA) for each observation with typical 
uncertainty 0.75 ms. TOAs were fitted to a timing model (see §3) using the pulsar timing 
software package TEMPO ^ 



3. Phase-Coherent Timing Analysis 

To determine spin parameters for PSR B0540— 69, we phase-connected 7.6 yr of timing 
data, that is, accounted for each turn of the pulsar. This is achieved by fitting TOAs to a 
Taylor expansion of pulse phase, 0. At time, t, (p can be expressed as: 

(f){t) = (Pito) + - to) + - + - tof + ■■; (2) 

where the subscript denotes a parameter evaluated at the reference epoch, to. TOAs and 
initial parameters are input into TEMPO, which gives as output refined spin parameters 
and residuals. To determine an approximate input ephemeris, we performed a periodogram 
analysis on 5 observations from MJD 51197 - 51206. This provided 5 values for u, on which 
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we performed a least-squares fit, resulting in an initial estimate for v. We then applied the 
center value of v and the fit value of v to the phase data and 'bootstrapped' a more accurate 
ephemeris which we extended over long time spans. 

We performed two timing analyses on the data. In §3.1, we describe our partially 
phase-coherent analysis, which determined spin parameters while minimizing the effects of 
timing noise and covariances arising from a polynomial fit to the data. In §3.2 we describe 
our fully phase- coherent timing analysis, done to confirm our partially coherent results and 
examine the timing noise apparent in the data. §3.3 describes a position fit performed using 
phase-coherent timing. 

3.1. Pcirtially phcise- coherent timing analysis 

The traditional method of phase-connecting all available data has resulted in many im- 
portant pulsar timing results owing to the accuracy available with absolute pulse numbering. 
However, in cases where timing noise seriously contaminates the data, many higher-order 
frequency derivatives are required to satisfactorily fit the data by removing red noise. This 
invariably introduces serious covariances between the fit spin parameters, without entirely 
removing the contamination from the timing noise. 

A partially phase- coherent analysis is less sensitive to contamination from timing noise. 
In this method, we phase-connect subsets of the data, where the length of each subset is 
determined by it being the longest that requires only v and i> be fit while having residuals 
consistent with white noise. In this way, the 7.6 yr of data were divided into 22 subsets. The 
measured values of v in these subsets are shown in Figure 1. A discontinuity can clearly 
be seen in this plot. We interpret it as a ghtch that occured at MJD 51342±24, where 
the epoch is taken to be the midpoint between the two neighbouring measurements, and the 
uncertainty is the range between these values. In §3.2 we confirm this interpretation using the 
more traditional fully coherent fit. A more accurate glitch epoch is impossible to ascertain 
from the partially coherent method due to sparse sampling and the small magnitude of the 
glitch. 

A value for v obtained from the slope of the best-fit line to all of the data will clearly 
be different from that obtained in two separate fits of the pre- and post-glitch data. The 
pre-glitch value of v implies n = 2.135 ± 0.016, while the post-glitch value yields n = 
2.144 ± 0.007, where the uncertainties are obtained by a bootstrap analysis which gives 
accurate uncertainties in cases where the formal uncertainties are thought to underestimate 
the true values, i.e. in the presence of timing noise (Efron 1979). The two values of n are 
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in agreement within the quoted uncertainties; thus we quote the average of the pre- and 
post-ghtch values, n = 2.140 ± 0.009. 

Figure 2 shows the measurements of z> with the pre-ghtch slope fitted out. This plot 
clearly shows the change in at the epoch of the glitch. The change in Az>/z> is ~ (1.5 ± 
0.1) X 10~^. The significant trend in the post-ghtch data in Figure 2 can be readily explained 
by timing noise. 

By further subdividing the data, we measured 38 values for i/. Fitting the pre-glitch 
trend (z> and i>) and removing it from the data shows the influence of the glitch on u as 
presented in Figure 3. No jump in i> can be seen at the glitch epoch owing to the small 
magnitude of the glitch (as shown in §3.2). In fact, Ai/ is on the same order as the uncertainty 
in the frequency measurements, rendering the partially phase-coherent method unsuitable 
for measuring a meaningful Au/u in this case. Spin and glitch parameters obtained in the 
partially coherent analysis are given in Table 1. 

3.2. Fully Phase-Coherent Timing Analysis 

To confirm the result obtained with the partially coherent timing analysis and obtain 
more precise glitch parameters, we phase connected all 7.6 yr of timing data. We first con- 
nected the data without allowing for a glitch in order to see if one was required by the fit. 
No phase jump was apparent at the epoch of the reported glitch, MJD 51325, however, we 
found that the timing residuals grew very significantly when connecting over this epoch. In 
order to phase connect the full data set with frequency derivatives and no glitch, u and 11 
frequency derivatives were required, as shown in Figure 4, resulting in n — 1.968 ± 0.004. 

We then allowed for a glitch near the epoch of that reported by Zhang et al. (2001). 
In this case, the data can be fit well with 3 glitch parameters (epoch, Au, and Au) and 3 
spin parameters {u, i>, and i>). Residuals for this fit are shown in the top panel of Figure 
5. The best fit glitch parameters, as determined with TEMPO, gives a ghtch epoch of MJD 
51335 ± 12, in agreement with the partially phase-coherent analysis. 

We found Au/u — (1.4 ± 0.2) x 10~^, which was undetectable in the partially coherent 
analysis, and 2.5(7 from the value reported by Zhang et al. (2001). We also found Az>/z> = 
(1.33 ± 0.02) X 10~^, in agreement with the partially coherent analysis, but lOcr larger than 
that reported by Zhang et al. (2001), Az>/z> = (0.85 ± 0.05) x 10"! 

To conclusively show that the model with the ghtch is a better description of the data, 
we fitted additional frequency derivatives until the same number of parameters was included 
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in each fit (i.e. 12 parameters). The bottom panel of Figure 5 shows the results of this 
fit. Note that the RMS residual is a factor of 5.15 smaller for the model with the glitch, 
corresponding to a reduced of 1.86 for 545 degrees of freedom, compared to a reduced 
of 49.4 for 545 degrees of freedom for the fit with no glitch. This confirms the conclusion 
of our partially coherent analysis, namely that a glitch occured as reported by Zhang et al. 
(2001). 

The phase-connected solution with the minimum number of derivatives (i.e. the number 
required to obtain phase-connection, in this case, 2), and the glitch fitted, gives a value of 
n = 2.10657 ± 0.00001. However, since there is clearly low-frequency noise contaminating 
the data (top panel of Figure 5), the formal uncertainty greatly underestimates the true 
uncertainty on n. Spin parameters vary when higher order derivatives are fitted due both 
to timing noise and covariances between parameters (Livingstone et al. 2005). Thus, we 
estimate the upper limit on the uncertainty in n, albeit not rigorously, from the variation 
in the measured value as higher order derivatives are added to the fit. This method implies 
an uncertainty of 2.11 ± 0.06, in agreement with the value obtained in the partially phase- 
coherent analysis. All spin and glitch parameters from the fully phase-coherent analysis are 
given in Table 1. Uncertainties given for parameters are the formal uncertainites, except in 
the case of n. 



3.3. Timing Position 

Our initial timing analysis was first performed using the pulsar position determined 
in a Chandra observation (Kaaret et al. 2001). The reported la uncertainty was 0.7" on 
the Chandra position. Holding this position fixed, the post-fit residuals for the fully phase- 
coherent analysis showed periodic behaviour, with a period of 1 year and amplitude ~ 
5 ms, shown in the top panel of Figure 6 (which also has fitted the glitch and 8 frequency 
derivatives as in Figure 5). This suggested a small but highly significant position hence 
barycentering error. We found that the periodicity in the residuals changes negligibly when 
higher order derivatives are fitted, implying that a polynomial does not describe the trend 
in the data. We performed a grid search for position by re-barycentering the data and 
calculating TOAs for each new position. We then phase-connected each set of new TOAs 
and calculated post-fit residuals. We found that the position that minimizes the RMS 
residuals is RA=05M0"^lin6 ± 0^04 and Dcc= -69°19'53'.'9 ± 0'.'2, 1'.'3 from the Chandra 
position. The phase- coherent RMS residuals fall by a factor of ~ 2 and the reduced is 
smaller by a factor of 4 when the fitted position is used. More importantly, the periodic 
trend in the data is removed by this fit, shown in the bottom panel of Figure 6. Remaining 
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systematics are likely due to timing noise. All quoted parameters for the phase-coherent 
and partially phase-coherent analysis are calculated using the new position. However, none 
of the spin parameters is significantly different from what is obtained using the Chandra 
position. To our knowledge, this is the first time phase-coherent X-ray timing has been used 
to try to improve a source position. If correct, we have effectively demonstrated superior 
source localization using the non-imaging, 1°-F0V PCA compared with that attainable with 
Chandra. However, it is possible that the source timing noise is merely mimicking a position 
error; continued timing will confirm or refute this result. 



4. Discussion 

4.1. The Braking Index 

By performing both a phase-coherent and a partially phase-coherent analysis on 7.6 yr 
of timing data for PSR B0540-69, we have found n = 2.140 ± 0.009. This value, stationary 
across the glitch epoch, is similar to that reported by Deeter et al. (1999), n — 2.080 ±0.003 
based on five years of GINGA data, and 1.7cr from that obtained by Cusumano et al. (2003), 
n = 2.125 ±0.001, based on 4.6yr of RXTEdata. Cusumano et al. (2003) however, reported 
no glitch and their uncertainty on n does not account for the effects of timing noise. The 
reason for the good agreement between our measured values is that despite the fact that 
they report no glitch, the value given for n is obtained from phase-coherent fits to the data 
before and after the glitch reported by Zhang et al. (2001), instead of a fit to all 4.6 yr as 
would be appropriate if no ghtch had occurred. We phase connect the same 4.6 yr of data 
they used without a glitch and find n — 1.73 ± 0.03, significantly different from the value 
obtained not considering the glitch, further implying the occurance of a glitch. 

Our value of n is significantly larger than that reported by Zhang et al. (2001), n — 
1.81 ± 0.07. We find that if we phase-connect only the same 300-day data set they used, we 
obtain a value n = 1.82 ±0.01, in agreement with their result. Thus their measurement was 
clearly contaminated by timing noise and/or glitch recovery, and affected by the relatively 
small time baseline used to measured n. 

The value oin — 2.140 ± 0.009 is significantly less than that for simple magnetic dipole 
radiation, n = 3. In fact, all measurements of n thus far are less than 3. Several explanations 
for this have been proposed. Loss of angular momentum owing to a particle wind (which 
young pulsars are well known to have) could account for this (Manchester & Peterson 1989). 
Another explanation is that the neutron star loses additional rotational energy by torquing 
a disk of supernova fallback material (Menou et al. 2001). Another explanation is a time- 
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varying magnetic field (Blandford & Romani 1988). The presence of a time- varying magnetic 
field can be determined by measuring a value of the third frequency derivative larger or 
smaller than that predicted by the spin-down model. The simple spin-down model given 
in Equation 1 assumes a constant value of which is only valid if the magnetic field, 
moment of inertia, and angle between the spin and magnetic axes are all constant. With 
present data, the third frequency derivative is not measurable for PSR B0540— 69 because of 
contamination from timing noise. It is possible that, as observations continue, this parameter 
will be measurable. For both the Crab pulsar and PSR B1509— 58, values of v have been 
measured and are consistent with the spin-down law (Lyne et al. 1993; Livingstone et al. 
2005). 



4.2. Glitches and Timing Noise 

We found that the best description of the 7.6 yr of RXTE data include a glitch at MJD 
51335 ± 12. Our partially phase-coherent analysis shows a clear discontinuity in z> (Figures 
1 and 2). The fully phase coherent analysis confirms this result. This glitch was previously 
reported by Zhang et al. (2001), at MJD 51325 ± 45 and later refuted by Cusumano et al. 
(2003), who claimed that the data could be described by fitting higher order derivatives, that 
is, timing noise. The detected glitch has a small magnitude in both Az//z/ and Az>/z> which 
Cusumano ct al. (2003) cite as evidence that no glitch occurred. However, fitting derivatives 
to the data shows that the residuals increase dramatically when fitting over the epoch of the 
reported glitch. This is exactly what is expected from a glitch, and is not expected from 
timing noise, which increases residuals gradually as observations are added. In addition, the 
partially coherent analysis shows a clear discontinuity in z> near (within uncertainties) the 
epoch of the glitch found by the phase-coherent analysis (Figure 5). Sparse observations of 
this source previous to the RXTE timing program most likely would not have been able to 
distinguish a small glitch from timing noise. Undetected glitches of the same magnitude as 
the one reported here are a likely explanation of the wide range of reported values of n. 

Although the magnitude of the ghtch is small, it is similar to the smallest glitches expe- 
rienced by the Crab pulsar, which is close in age, spin-down luminosity {E = —An'^Iuu), and 
dipole magnetic field {Bq ~ 3.3 x 10^^ (PP)^/^) to PSR B0540-69. In fact, PSR B0540-69 is 
often referred to as the "Crab Twin" due to these similarities. The Crab pulsar has glitched 
with parameters as small as Au/u ~ 2 x 10^^ and Au/u ~ 10^^ (Lyne et al. 1993). However, 
the Crab pulsar also experiences glitches on the order of Av jv ~ 10~^ with values of Ai> jv 
as large as 4 x 10~^ (Wang et al. 2001). None of these larger magnitude ghtches has been 
observed in PSR B0540-69. 
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The glitch activity parameter, Ag, is defined as the cumulative change in frequency from 
all glitches over the observation length, = E(Az//z/)/At (McKcnna & Lyne 1990). The 
Crab pulsar has been observed to glitch 14 times between 1969 and 2000 and has a glitch 
activity parameter A^ ~ 0.1 x 10~^yr~^ (e.g. Shemar & Lyne 1996). PSR J1119— 6127, aged 
~ 1600 yr with a measured braking index, n = 2.91 ± 0.05 has been observed to glitch once 
implying A^ ~ 0.02xl0-^yr-i (Camilo et al. 2000). However, the fact that PSR J1119-6127 
glitched during less than 2 yr of observations imphes that this pulsar likely ghtches often and 
may prove to have a glitch activity parameter on the order of that of the Crab pulsar, or 
larger. At the opposite end of the spectrum hes the ~ 1700 year old pulsar PSR B1509— 58, 
which has not glitched in 21.3 yr of observations (Livingstone et al. 2005). 

The single observed glitch in PSR B0540-69 between 1996 and 2004 implies A^ ~ 
0.002 X 10~^yr~^. Undetected glitches occuring in the period prior to RXTE observations 
would increase A^ and could be the source of discrepant values of n in the literature. However, 
it is also possible that the pulsar did not glitch during this time period, which would lower 
the overall value of Ag. In this case, the range of reported values of n would be due to 
timing noise and/or pulse counting errors. The sparse nature of observations (and in some 
cases data quality) between 1979 and 1996 precludes conclusions on this point. Although 
the estimate of A^ is based on a single glitch and is therefore extremely sensitive to any 
future ghtch activity, it is clear that PSR B0540— 69 glitches less than the Crab pulsar and 
somewhat more than PSR B1509— 58. Thus, there is a range of possible glitch activity 
amongst very young pulsars. 

The variety of glitch activity is suggestive of a range of internal temperatures among 
young pulsars, since glitch activity is possibly associated with neutron-star temperature 
(McKenna & Lyne 1990). In this theory, the very young pulsars such as the Crab, PSR J1119— 6127, 
PSR B1509— 58 and PSR B0540— 69 have small ghtch activities because they are relatively 
hot. Pulsars that are slightly older, e.g. the 10-kyr old Vela pulsar, are cooler and have 
larger glitch activities, ~ 8 x 10~^yr~^ (Dodson et al. 2002, and references therein). Finally, 
as pulsars become much older and cooler, they are believed to eventually stop glitching en- 
tirely. However, recent reports of glitches in AXPs, which are thought to be hot neutron 
stars, argues that temperature is not the only factor in determining glitch behaviour (Kaspi 
et al. 2000; Kaspi & Gavriil 2003; Dall'Osso et al. 2003). Regardless, it is intriguing that 
there is such a wide variety of glitch behavior in young pulsars. 

This research made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center Online Service, provided by the NASA-Goddard Space Flight Cen- 
ter. VMK is a Canada Research Chair and an NSERC Steacie Fellow. Funding for this 
work was provided by NSERC Discovery Grant Rgpin 228738-03 and Steacie Supplement 
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Table 1. Parameters for PSR B0540-69. 
Parameters for timing analysis 

Dates (Modified Julian Day) 50150 - 52935 
Right Ascension^ (J2000) 05^40"'11U6 ± 0!04 

Declination^ (J2000) -69°19'53'.'9 ± 0'/2 

Parameters for partially coherent analysis 

Braking Index, n 2.140(9) 
Glitch Epoch (MJD) 51342(24) 

Au/u (10-^) 1.5(1) 

Parameters for phase-coherent analysis 

Epoch (Modified Juhan Day) 51197.0 

z/ (Hz)*^ 19.80244383176(2) 

z> (10-1° s-2)b -1.878039597(8) 

i) (10-21 s-3)b 3.752027(2) 

Braking Index, 2.11(6) 

Glitch Epoch (MJD) 51335(12) 

Au/u (10-9)^ 1.4(2) 

A;>/;> (10-7' 1.33(2) 



^Position fit from a grid search by minimizing the 
phase residuals in a phase-coherent analysis, see §3.3. 

^Quoted uncertainties are formal uncertainties as 
reported by TEMPO. 

^Uncertainty determined from variation in n as 
higher order derivatives are fitted, as explained in 
§3.2. 
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Fig. 1. — Measurements of z>, the slope is i>. The ghtch occuring near MJD 51342 is shown 
with an arrow. The pre-ghtch value isi> — (3.81±0.03) x 10~^^ implying n — 2.315±0.016, 
while the post-glitch value is i) = (3.81 ± 0.01) x 10"^^ s'^ implying n = 2.144 ± 0.007. The 
average of pre- and post-glitch n is 2.140 ± 0.009. Measurement uncertainties are smaller 
than the points and are omitted. 
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Fig. 2. — Measurements of v with the pre-gUtch slope (ie. v) subtracted. The value of 
IS.V jv ~ (1.5 ± 0.1) X 10~^ is in agreement with the value obtained from a phase-coherent 
analysis. The vertical line indicates the approximate ghtch epoch, MJD 51342. The trend 
in the post-glitch data is attributable to timing noise. 
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Fig. 3. — Measurements of u, with the pre-ghtch trend removed. The arrow marks the ghtch 
epoch, MJD 51342. Uncertainties are smaller than the points. 
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Fig. 4. — Residuals for PSR B0540— 69 with v and 11 frequency derivatives - and no glitch 
- fitted. This was the minimum number of derivatives needed to connect the data, i.e. yield 
residuals having amplitude less than 0.5 P. 
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Fig. 5. — Timing residuals for PSR B0540— 69. Top panel shows residuals with three glitch 
parameters and two frequency derivatives fitted. This fit yields a value of n = 2.11 ± 0.06. 
Bottom panel shows residuals with 12 parameters total fitted, for comparison with Figure 4. 
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Fig. 6. — RMS residuals for the Chandra position (top) and the best-fit position of 
PSR B0540— 69 (bottom). Both sets of residuals have had u and 8 frequency derivatives 
fitted; removing higher order derivatives does not improve either fit significantly. 



